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ABSTRACT 

o 

Aims. An analytical solution for the discrepancy between observed core-like profiles and predicted cusp profiles in dark 
■ matter halos is studied. 

Methods. We calculate the distribution function for Navarro-Frenk- White halos and extract energy from the distribu- 
^"^ , tion, taking into account the effects of baryonic physics processes. 

Results. We show with a simple argument that we can reproduce the evolution of a cusp to a flat density profile by a 
decrease of the initial potential energy. 



o 

C"| I There are many dwarf galaxies in the Universe (e.g. iMarzke fc da Costa (11997): ISawala et al.l (l2010f)h and the Local 
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1. Introduction 



Q^i Group is an excellent place to study these systems and their properties (|Mateolll998t IPasetto et al.H2010f ). 

The n early solid body rotation curve observed in dwarf galaxies indic at es a central co re in the dark matter distribu- 
2 ! tion (e.g.. iBurkertl (|l995[ k Ide Blok fc Bosrnal ([2002b : iGentile et al.l (|2005lh Ide Biol (12005D). These results disagree with 

■ predic t ions from numerical s imulations, which require density profiles with cusps (e.g. INavarro "eTaD(fl997l) :lM oore et al.l 
^ ; (fl999l) : INavarro et al.l dgOOl ) . 

— i! INavarro et al.l (|l99ot) were one of the first to point out the possibilit y that feedback mecha nisms can turn the central 

■ dark- matter cusp into a cored one. Following this line of investigation, iGnedin fc Zhao! (|2002l ) analyzed the influence of 
winds. iRead fc Gilmord (120051) discussed how extern al impulsive mass- loss events can successfully act as a flattener for the 



■ central density cusps, and Mashchenko et afl ([20061 ) showed that random bulk motions of gas in small primordial galaxies 
' could flatten the central dark matter cusp in ~ 10 8 years. Beyond that, several authors have suggested that the interstellar 
medium (ISM) of dwarf galaxies systems c ould be entirely removed as a result of sta r formation (iDekel fe Silk 1986; 
Mori et all 119971: iMurakami fc Babul Il999t iMac Low fc Ferraral [1991 iMori et al.l [20021: iHensler et all 12004 iMori et al l 



20041) . 



m 
in 

In the present work, we construct a simple model under the hypothesis that astrophysical processes are able to 
f— ^ remove baryonic gas from dwarf galaxies. We show for the first time a clear and analytical argument to solve the cusp 
I | problem. We argue that baryonic feedback could change the distribution function on the phase space of energy and 
• • angular momentum, flattening the density profile. 

This paper is organized as follows. In Sect. 2, we briefly discuss the removal of baryonic gas from dark matter halos 
, as a consequence of different astrophysical processes. In Sect. 3, we review the formalism of distribution functions. In 
; i ■ Sect. 4, we discuss the changes in density profile caused by gas removal and show our results. Finally we present our 
' conclusions in Sect. 5. 



2. Expulsion of baryonic gas 

Our main goal is to give an analytical explanation to the flat density profile observed in dwarf galaxies. Our starting 
point will be the results obtained from ACDM numerical simulations ([Moore et a.1.1 119991: INavarro et al.l I2004T) . In this 
context, we begin with a dwarf galaxy that has a Navarro-Frenk- White profile (NFW) (| Navarro et al.|[l997f) after the 
formation of its first stellar population, and study the signatures left in the dark matter profile when part of its baryons 
are removed. The first natural question is: why are the baryons removed? 

Before answering that question, it is important to emphasize that a cusp profile is not only characteristic, but 
necessary. In the initial phase of galaxy formation, when only a sm all part of the stellar population is forme d, a deep 
gravitational potential well is needed in order to trap the baryons ([Mori et al.l 119971 : iBradamante et aT1ll998l) . Beyond 
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that, baryon accr etion in the center of the halo will suffer adiabatic contraction and drag the dark matter profile into an 
even steeper one (Blu menthal et al.lll986h . 

On the other hand, the baryonic characteristic of the gas causes it to radi ate energy, increasing the entropy of the 



dark matter background and leaving strong signatures in the density profile (jEl-Zant et al.l 12001 : Romano- Diaz et al.l 



2008). In this context, a lso the effects of gas bulk motions need to be considered, as shown in ' iMo fc Mad (|2002h and 



Mashchcn ko et al.l (2006) 



Af ter th is initial collap s e, the galaxy will go through one (or probably more) star-formation bursts . [Mac Low fc Ferraral 
(1999) and iPasetto et al.l (l2010D showed that SF alone c an rem ove almost all the interstellar baryonic gas inside dwarf 



galaxies, for M^aio < 10 7 M ? . iGelato fc Sommer-Larsenl (|1999D explored how different m ass ejection fra c tions and disk 



density distribution could fit the observed rotation curve of the dwarf galaxy DDO 154. lHavashi et al.l (|2003l) demon- 
strated that the mass loss can vary from 5% to 95% as a result of tidal effects caused by a host larger galaxy. 

Below we will consider only astrophysical effects that are able to remove fr om a small percentage up to all the baryonic 
gas initially trapped by the dark matter halo. Following IPasetto et al.l (|2010f) . we estimate that 1/50 of the total mass in 
dwarf galaxies is visible. Moreover, 15% of this visible mass is in form of stars and 85% is in form of interstellar gas. We 
consider that the astrophysical processes cited before are able to remove from 5% to 100% of the gas in the galaxy halo. 
This means that for halos with M\ la i ~ 10 7 M Q , from 0.085% to 1.7% of the total halo mass can be removed. 

Note that the crucial point in our analysis is the percentage of total halo mass lost. All the other astrophysical 
mechanism cited above that do not necessarily remove baryons are not considered here. Our results are valid for any 
mass loss mechanisms, as long as they expel the same percentage of baryons out of the dark matter halo potential well. 

3. Dark matter density profile 

At a given time, t, a full description of any collisionless system can be obtained by specifying the number of particles, 
f(x,v,t)d 3 xd 3 v, that have positions in the small volume d 3 x centered in x and small r ange of velocities d 3 v centered in 
v. f(x,v,t) is called the distribution function, hereafter DF (jBinnev fc Tre mainc 200^). The DF of a mass distribution 
in a steady state is related to its density profile through 



p(r) = I f(r,v)d 3 v . (1) 

For any general spherical system, we can defi ne a DF, f(E, L), which depends on the energy, E, and the angular 
momentum absolute value, L. lOsipk ov (1979]) and iMerrittl (|1985|) found that in anisotropic spherical systems, / can be 
expressed as: 

f(Q,L 2 )=T(Q)L 2a , (2) 



where a is a real number and 



for a radial parameter r a , known as the anisotropy radius. It is useful to define the relative potential, 'J = — 0, where <f> 
is the gravitational potential. 

For distribution f unctions that can be writt en in the form of Eq. ([2]), the fundamental integral equation can be written 
as (lCuddefordlll99ll:lBinnev fc Tremaindl2008l) 



Ztt 



r 



p(r) = — / F{Q) dQ (4) 



x 



2r 2 (*-Q)/(l+^-) L 2a dL 2 



/2(*-Q)-(£Vr- 2 )(l + f) 
After performing the integration over L, we obtain a single integral equation for J-(Q), 

p(9) = Q(a) f (* - Q) a+1/2 F{Q)dQ , (5) 
Jo 



where 



6(a) = 2 "+ 3 /2 7r 3/2 r(a + l) 2a + _ (6) 



r(a + 3/2) V t. 

For a general a, Eq. ([5]) can be inverted, resulting in 



7re(a)(a + l/2)!dQ»+iy (Q _*)«+!-» ' K } 
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with n defined by 

ra=[a + l/2] + l, (8) 

and [a + 1/2] representing the largest integer less than or equal to (a + 1/2). 

From the above expressions, it is clear that any perturbation in E and/or L will change the DF (Eq. ©) and, 
consequently, leave a signature in the galaxy density profile (Eq. (JTJ). Our strategy in this work is to start with a NFW 
profile in Eq. (J7J, find it a suitable DF and analyze the change in the density profile caused by perturbations in the DF. 
Although, as we initially proposed an entirely analytical calculation, we must express the NFW profile in a way that 
makes a completely analytical approach possible. This is done in the next section. 

4. Erasing the cusp 

4.1. Spherical density profiles 

In order to evaluate the distribution function for a cusp-like profile, we need to consider an initial densit y profile with a 
centr al cusp and place it into Eq. (0 . For simplicity, we adopt the following family of spherical potentials (jRindler-Dallerl 



Si, 



V Kb/ J 

where <j)Q is the value of potential inside the halo radius, 77, 7, and b are positive parameters. Using the Poisson equation, 
we found the dimensionless density profile corresponding to this potential, 

~ { ,_ b^+ 2 (1 + 77)6" + (1 - 777)^ 

P[T> 1 + ry r 2 -v(bv + r")^+ 2 ' 1 ' 

Defining = — </>(r)/</3o, we can also write the density as function of the relative potential, 



p(^f) = — ( [(1 - 777)* * +1 + 77(7 + 1)*" 
T] + 1 I- L 



^+1 



x 



*-v7_i) n, (ii) 



where the tilde denotes dimensionless quantities (for more details see iRindler-Dalierl (|2009h ). The Eq. (HQ) is general 
enough to be applied to any spherical density profile, given an appropriate choice of parameters. The advantage of this 
expression is that it permit us to write p in terms of ^P, which is not always analytically possible. However, we can always 
find a fit of Eq. (fTTj) corresponding to a specific density profile. 



4.2. Connection to the Navarro-Frenk-White profile 

One o f the strongest pred ictions from cosmological simulations is the universal density profile present in dark matter 
halos (jNavarro et al.lll997l ). 

p(r) = . f f s . .„ , (12) 
KV ' r/rs{\ + r/rs) 2 v ; 

where p s is de central density and r s is the scale radius of the halo. 

Using 77 = 1, b = 1 and 7 = 1/2, Eq. (ITU1) becomes numerically similar to NFW profile, as we can see in Fig. [T] Given 

this choice of parameters, \&, Eqs. (fTU|) and (jTTJ) become 

*(r) = -7^=, (13) 
vr + 1 

^ 4,(r + lV (U) 
and _ _ 

3^7,^,5 

= a — aTT' ^ 

4 - 4^ 2 

In order to test our calculations, we first found a DF for the NFW profile and reconstructed the density profile using 
the Eq. ([13) in Eq. ©. We found that for a - 1/2 and r a - 1, Eq. © agrees well with the DF for NFW profile. With 
these parameter values we have 

_ _Q>(UQ> -3W^t«| . ( 

2tt 2 (Q 2 -1) 3 
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0.001 0.01 0.1 1 

Fig. 1. Reconstructed density profile (red dashed line) compared with original NFW profile (blue solid line). 

The above expression placed into Eq. §5§ recovers Eq. (jT5j) . 

We suppose that the macroscopic variables, like radial density profile, will change if the halo experiments a physical 
process that perturbs its DF. We considered that the removal of baryonic gas will shift the energy of the system by 
Q — > Q ± AQ. According to the arguments in Sect. 2, SF and external tidal effects can remove (0.085 — 1.7)% of the halo 
gravitational potential energy. In order to find a new DF configuration, we have to evaluate J-(Q — AQ) « J-{Q) — AQ§f^. 

Which is just a solution of the Boltzmann equation, ^ + = 0, if we consider that the gas removal occurs in a 

short time interval, At ss dt. 

The results are presented in Fig. [2j where we show the new configuration of the density profile after remotion of 
(0.085 — 1.7)% from gravitational potential energy. It is evident from this figure that a little change in the total energy 
of halo can flatter the profile in the inner region. Note that the core-like behavior we were looking for is already evident 
for AQ 0.001. As a consequence, any feedback capable of removing ss 10% of interstellar gas would be able to smooth 
the matter density profile. 

5. Conclusions 

The N-body simulations based on collisionless cold dark matter do not show the core-like behavior observed in local 
dwarf galaxies. Instead, they are better described by a steep power-law mass density distribution called cusp. Several 
attempts were made with numerical simulations in order to take into account the gas dynamics in the evolution of dark 
matter profiles. 

Considering that the total energy of the halos will change because of baryonic feedback, we showed that it is possible 
to model the evolution of an initial cusp-like profile into a core-like one. We used an analytical analysis based on the 
density profile DF. 

The advantage of this approach is that we obtain resu l ts similar to those extra cted from numerical simulations, which 
take into account the gas dynamics (|Pasetto et al.|[20Tot iGovernato et al.|[2010D in a very simple way. On top of that, 
considering the baryon loss as a mechanism of flattening the dark matter density profile of dw arf galaxies allows us to 
support the low amount of baryons observed in dwarf spheroidal galaxies ( Gil more et al.l[2b07t ). 

However, as discussed in the introduction, there are two competitive processes occurring during the formation of 
galaxies, baryon accretion in the center of the halo will suffer adiabatic contraction and drag the dark matter profile into 
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Fig. 2. New equilibrium configuration of the perturbed halo removing 0.1% (red dot-dashed line) and 1.7% (dashed green 
line) of total energy Q, compared to the NFW profile (blue solid line). 

an even steeper one, and on the other hand, the baryonic characteristic of the gas causes it to radiate energy, increasing 
the entropy of the dark matter background and preventing the contraction. As our knowledge about these two processes 
is still incomplete, it is important to take into account these effects in future works. 
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